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Abstract 

A popular derivation of the apparent temperature of the particle transverse mo¬ 
mentum spectra and its inconsistency are considered. An improved formula for the 
apparent temperature is discussed. It is shown that recent results on transverse mass 
spectra of J/ip and ip' mesons support a hypothesis of statistical production of charmonia 
at hadronization and suggest the early thermal freeze-out of these mesons. Using the 
apparent temperature formula the collective transverse velocity of the hadronizing quark 
gluon plasma is extracted to be (vp) ~ 0.2. Predictions for transverse mass spectra of 
hidden and open charm mesons along with bottomonium at SPS and RHIC energies are 
discussed. 

1. Statistical Production of Particles. J/ip suppression in nuclear collisions was sug¬ 
gested by T. Matsui and H. Satz to be a signal of the deconfinement phase transition. Nowa¬ 
days there are several models which mainly differ from the original idea of the J/ip suppression 
in the mechanism of charmonia dissociation either in quark gluon or in nuclear media. Re¬ 
cently a principally different picture of the statistical charm production was suggested JTJ . The 
next breakthrough in this field, the statistical coalescence model, suggests J|] that charmo- 
nium is generated by the coalescence of earlier produced c and c quarks. However, important 
questions, what happens after hadronization and what are the transversal mass spectra of 
J/ip and ip' mesons, were not considered in J^]. Therefore, I will mainly concentrate on 
the problem of the J/ip and ip' freeze-out (FO) and will analyze their apparent temperature 
according to Q. Then I will answer the question, “What can we learn from experimental 
data about hadronization of quark gluon plasma (QGP)?” 

The idealized concept of chemical (hadron multiplicities) and thermal (hadron momentum 
distributions) freeze-outs enables one to interpret data on hadron production in relativistic 
nucleus-nucleus (A+A) collisions. The first experimental results on yields and transverse 
mass (rriT = yjpp + m2 ) spectra suggest the following scenario: for the most abundant hadron 
species ( 7 r, N , K , A) the chemical FO , which seems to coincide with the hadronization of the 
QGP, is followed by the thermal or kinetic FO occurring at a rather late stage of the A+A 
reaction. Thus for the central Pb+Pb collisions at 158 A-GeV the temperature of the chemical 
FOwas extracted from the fit of the multiplicity data to be Tp = 175 ± 10 MeV |]J. The 
thermal FO parameters, i.e., temperature Tp and averaged transversal velocity (vp ), for pions 
were determined from the results of two pion correlations || to be quite different 

T f = 120 ± 12 MeV, (v T ) = 0.5 ± 0.12 . (1) 

On the other hand from hadronic cascade simulations it is known that thermal FO of 
multistrange hadrons (0, S, Q) happens, probably, earlier than the kinetic FO of pions. Now 
the main question is, “When does the J/ip freeze-out occur? Late or early?” 

The convolution of the transverse flow velocity of the matter element with the thermal 
motion of hadrons in the rest frame of this element leads to a nearly exponential shape of 
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final rriT spectrum ^ ~ C ■ e m T/T with the apparent temperature (AT) T* which 

is defined from the distribution at hxed longitudinal rapidity y p of the particle as 
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In the limit mx » T and px 


dpf \ dUpdpf, 

0 the AT depends on the particle mass and (vx) as follows 

T* = Tp + a ■ m ■ ( vx ) 2 • (3) 

For the spherical hreball Eq. (|3|) was hrst found in Ref. For the cylindrical geometry, 
surprisingly, there are two answers which differ by the value of the coefficient a: there is a 
naive result a — \ 0 and more elaborate one ot—\ [§]• Therefore, before going further it is 
necessary to discuss the origin of the difference between the results of Refs. 0 and ||. 

2. Particle Spectra and Apparent Temperature. The general expression for the 
momentum distribution of the outgoing particle having the 4-vector of momentum p p = 
= (rrir cosh y p , p T cos( 0 p ), pxsm((f p ), rrix sinh y p ,) which is emitted from the fluid element 
having hydrodynamical 4-velocity u p = 7 x(cosh?/L, vt cos(</> n ), Vt sin(0 U ), sinh yL ,) of the 
arbitrary FO hypersurface £ is given by the cut-off distribution function || [T0| 
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where the standard notations are assumed for 4-vectors (y p , yx and yL are particle longi¬ 
tudinal, fluid transversal and fluid longitudinal rapidities, respectively, and f denotes the 
corresponding polar angles), yr = 1 /\J 1 — vf is the relativistic 7 -factor, denotes the 

4-vector of external normal to the FO hypersurface £ and is the one-particle phase space 
distribution function. The O-function in (£|) is of crucial importance because it makes sure 
that only outgoing particles are counted from both the time- and space-like FO hypersurfaces 
which are defined, respectively, by the positive sign ds 2 > 0 and negative sign ds 2 < 0 of the 
element square ds 2 = dt 2 — dR 2 — dz 2 of the hypersurface (see 


TOU f° r detail). 


For further analysis I shall neglect the contributions coming from the time-like parts of 
the FO hypersurface and then comment on general case. Now the product p p d£ p > 0 is 
positive for any momentum and the cut-off distribution (f|) automatically reproduces the 
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famous Cooper-Frye result 0- For the Boltzmann distribution <p = 
emitted from the FO hypersurface R* = R*(t,z ) one obtains from (|]) 
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Now it is clearly seen that the main contribution to the integral (|6]) corresponds to the small 
angles 4> p between a 3-momentum of particle and a 3-vector of the hydrodynamic velocity, 
whereas contributions coming from the large angles f p are suppressed exponentially. 

Neglecting the term with in Eq. ([]) for px « mx and expanding the Bessel function 
Iq{x) « 1 + (0.5a:) 2 for pxVx7 t << T, one Ends from (^) the following result for the AT 
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where the double averaging means a double integration over time and longitudinal coordinate 
with the weight function defined by Eq. (^). For the longitudinal expansion depending on 
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the rapidity only (e.g. Bjorken expansion) the integration over ^-coordinate can be done 
because for heavy particles the Boltzmann exponential behaves like a Kronecker 5-function, 
i.e., e - m T7T cosh( y p - y L )/ T ^ This leads to a single averaging over the evolution time 

(last equality in Eq. ([?]) ) with a slightly modified weight function. The next approximation 
in ( 0 ) implies the nonrelativistic transversal expansion and m-r —* m 
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where the last equality is fulfilled within 10 % for the linear dependence of the velocity on 
the transversal radius. Note, however, that a popular expression (^) was obtained under 
the condition which is hardly fulfilled in practice for heavy particles and, hence, Eq. (||) 
cannot be established from Eq. (0)- Moreover, for particles of rrir > 1 GeV, freezing out 
under conditions (0), Eq. (0) leads to negative contributions to the AT and even to negative 
values of the AT! The latter follows from the fact that contributions coming from the small 
transversal radii (and, hence, small vt) in ( 0 ) are suppressed. 

Thus an accurate examination of the derivation of Ref. [|7] shows that negative AT can 
be seen, if the FO hypersurface is a space-like one. In Ref. [dQ] there are examples given 
that negative AT exist, if the Cooper-Frye formula is applied to the time-like hypersurfaces, 
whereas the cut-off formula (j4j) generates positive AT in that case. Discussion of the AT for 
the arbitrary hypersurface is out of the scope of this talk. I mention only that negative AT 
not seen experimentally may be an artifact of the Cooper-Frye formula. 

In contrast an improved derivation of Ref. [Q accounts for the finiteness of the system and 
is partly free of the problems discussed above. In j0] the Cooper-Frye formula is modified by 
an additional factor exp{—a(cosh?/ T — 1 )} (with a = 1 /((R*) u^( 0 )) 2 ) which in case of non¬ 
relativistic transversal expansion reduces to exp{— 2 h R ~. : i }• Here ( R*) is the mean transverse 
radius and i s Hie transverse velocity derivative at the center of the fireball. The large 

radii contributions in (|7|) are suppressed due to such a factor and, hence, they do not lead 
to negative AT values. Further evaluation ends in (||) with a = - which appears from the 
additional factor and reflects the cylindrical symmetry. Therefore, the latter is used below. 

3. Freeze-out of J/ijj Meson. Recently it was found that the data on J/'ij) and ip' yields 
in central Pb+Pb collisions at 158 A-GeV are consistent with the results of the statistical 
model for the typical value of Th = 175 MeV extracted from light hadron systematics (see 
Refs, in 10 )- The hypothesis of statistical J /ip production at hadronization can be further 
tested using data on m T spectra. New NA50 data [HJ on transverse mass spectra of J/ip 
mesons in central Pb+Pb collisions at 158 A-GeV confirm this expectation: the spectrum is 
nearly exponential with T*(J/ip) = 245 + 5 MeV. The measured T*(J/ip) value is significantly 
smaller than expected one from Eq. ([!) for the pion FO parameters (Tjj) (T* = 610 MeV). 

The ‘low’ value of T* for J/ip suggests its rather early thermal FO . The low interaction 
cross section of the J/ip meson with the most abundant hadrons [|E| and its large mass lead 
to a very low probability of the J/ijj rescattering on hadrons. Following 0, it is assumed 
therefore that the thermal FOof J/'ip coincides with the hadronization of QGP, i.e., that 
the J/ip meson does not participate in the hadronic rescattering after hadronization. It is, 
however, natural to expect that there is a significant collective transverse flow of hadronizing 
QGP developed at the early stage of partonic rescattering. Consequently the AT of the J/ip 
meson as well as all other hadrons for which chemical and thermal freeze-outs coincide with 
hadronization can be expressed as: 
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where (v^) is the mean transverse flow velocity of the QGP at the hadronization. Assuming 
Th = 175 MeV and using the measured value of T*(J/0) = 245 MeV one finds from Eq.(||): 
(vt) — 0.19. As expected the obtained transverse flow velocity of QGP at hadronization 
is significantly smaller than the transverse flow velocity of pions (~ 0.5). The linear m— 
dependence of (|9j) is shown in Fig. 1 by the lower solid line. Within the approach discussed 
here, Eq. (^) can be used to obtain a next estimate of the lower limit of the measured AT 
for all hadrons. In fact the values of the parameter T* for all light hadrons are higher than 
The recent results [19) on the my spectra of the 0' meson indicate that T*(0') & 
ThW) = 258 ± 5 MeV, which suggests that the 0' meson (like J/0) does not participate in 
the hadronic rescattering either. 

One may expect that the thermal FO may coincide with hadronization also for D and T 
mesons. Under this assumption one can calculate the value of the apparent temperature for 
those mesons : T*(D) = T^(D) = 217 MeV and T*( T) = T#(T) = 392 MeV, respectively. 

For the production of open and hidden charm particles in A+A collisions at RHIC energies 
one can expect stronger transverse collective flow effects than at the SPS. This will lead to a 
linear mass dependence (|]) of the AT with approximately the same value of T H = 175 MeV, 
but with a larger value of (v^). A recent analysis [^(J of the RHIC data leads to the above 
value Th = 175 MeV of the chemical FO temperature. 

Since the RHIC data for J/0 mesons are not measured yet, we use the hydrodynamic 
calculations of Ref. [2l]] which predict the value of (wff) = 0.30 at the hadronization in 
Au+Au collisions at RHIC. This leads to an increase of the AT of charmed hadrons at 
RHIC in comparison to those values at SPS, e.g., T*(J/0) = 20)(J/0) = 350 MeV and 
T*(D) = Tjj{D) = 280 MeV. If bottomonium experiences the same FO conditions at RHIC 
like J/0 meson, then its apparent temperature will be about T*( T) = T) = 720 MeV. 

In summary, recent results on transverse mass spectra of J/0 and 0' mesons in central 
Pb+Pb collisions at 158 A-GeV are considered. It is shown that these data support the 
hypothesis of the statistical production of charmonia at hadronization and suggest a simulta¬ 
neous hadronization and the thermal FOfor J/0 and 0' mesons. Based on this approach the 
collective transverse velocity of hadronizing quark gluon plasma is estimated to be (vj) ~ 0 . 2 . 
Predictions for transverse mass spectra of hidden and open charm along with hidden bottom 
mesons at SPS and RHIC are discussed. 
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and (uy) = 0.19 which correspond to the QGP hadronization. 
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